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In the region where the Higgs energy is large the process e + e _ — > ttH is governed by nonrelativistic QCD dynamics 
and one has to employ effective theory methods to make first principles QCD predictions. In this talk we use the 
effective theory vNRQCD to compute the Higgs energy distribution at next-to-leading logarithmic approximation. It 
is shown that the corrections are particularly important for smaller cm. energies. 

o : 
o : 

^ 1. Introduction 

c : 

It is one of the major tasks of future collider experiments to unravel details of the mechanism of electroweak 
symmetry breaking (EWSB). In the Standard Model (SM) EWSB is achieved by the Higgs mechanism. The particle 
masses are generated by the Higgs field vacuum expectation value V = (V^Gf) 1 ^ 2 ~ 246 GeV arising through the 
t-H Higgs self interactions. The mechanism also predicts that Higgs bosons can be produced in collider experiments. 
While a Higgs boson with a mass smaller than 1 TeV can be found at the LHC, precise and model-independent 
measurements of its quantum numbers and couplings can be gained from the e + e~ Linear Collider. A crucial 
prediction of the Higgs mechanism is that the Higgs Yukawa coupling to quarks X q is related to the quark masses 
by m q — X g V. At the e + e~ Linear Collider the top quark Yukawa coupling can be measured from top quark pair 
, production associated with a Higgs boson, e + e~ — > ttH, since the process is dominated by the amplitudes describing 
Higgs radiation off the tt pair. This process is particularly suited for a light Higgs boson since the cross section can 
• then reach the 1-2 fb level. Assuming an experimental precision at the percent level, QCD and electroweak radiative 
corrections need to be accounted for in the theoretical predictions. In the approximation that top quarks and the 
Higgs are stable particles the Born cross section was already determined some time ago in Refs. [lj. For the 0(a s ) 
QCD one-loop corrections a number of references in various approximations exist On the other hand, the 

full set of one-loop electroweak corrections was obtained in Refs. 0, Q and also in Refs. 0. In Ref. @ a detailed 
• i-H , analysis of various differential distributions of the cross section cr(e + e _ — > ttH) can be found. 

A particularly interesting kinematical phase space region is where the energy of the Higgs boson is large and close 
to its kinematic endpoint. The tt pair then becomes collinear and flies opposite to the Higgs direction to balance 
the large Higgs momentum, see Fig. ^ For large Eh, on the other hand, the tt invariant mass Q 2 approaches 
4m(, Eh = 2^75 ( s + m H ~ Q 2 )j an d the top quark pair is nonrelativistic in its own center-of-mass (cm.) system. 
Because the Higgs is very narrow for a mass below the W + W~ threshold, strong interactions between the tt pair 
and the hadronic Higgs final state can be neglected. Thus, close to the Higgs energy endpoint the tt QCD dynamics 
is exclusively governed by the nonrelativistic physics known from the process e + e~ — * tt in the tt threshold region 
at y/s 2m t . In this regime the so-called Coulomb singularities cx (a s /v) n , with v = (1 — Am 2 / Q 2 ) 1 / 2 being the top 
quark relative velocity in the tt cm. frame, arise and require predictions using an expansion in a s and v rather than 
just a perturbative computation in the number of loops. 

This singularity structure is most easily visible in the Higgs energy distribution, da(e + e~ — > ttH)/dEn- While the 
Born distribution approaches zero for Eh — > E^ ax , da/dEn ~ v y|, the 0(a s ) fixed-order perturbative corrections 
are proportional to a s at the endpoint 0,0] and the 0(pt 2 s ) corrections even diverge like a 2 ./ v. The problem might 
be avoided by imposing a cut on Eh or Q 2 , but this is unnecessary because there exists an elaborate technology 
being used for the threshold region in the process e + e~ — > tt 8] that allows for systematic QCD predictions with 
renormalization group (RG) improvement. Imposing a cut would be also disadvantageous as the nonrelativistic 
portion of the ttH phase space increases for smaller cm. energies which are relevant for a measurement in the first 
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Figure 1: Typical constellation of momenta for the process e + e~ — > ttH in the large Higgs energy endpoint region. 



phase of the ILC program. Because the SM top width is quite large, T t w 1.5 GeV, the corresponding QCD effective 
theory computations can be carried out with perturbativc methods for all Higgs energies in the endpoint region. 

In this talk we present the Higgs energy distribution da/dEn in the large Higgs energy endpoint region at NLL order 
in the nonrelativistic expansion using the framework of "velocity" NRQCD (vNRQCD). For details on the conceptual 
aspects, concerning powercounting, the operator structure of the effective theory action, and renormalization we refer 
to Refs. For a more detailed discussion of the computations for this work see Ref. 0. 



2. Review of Effective Theory Ingredients 



The effective theory vNRQCD provides a systematic RG-improved description of dynamics of nonrelativistic tt 
pairs. The system is characterized, for any energy in the threshold region, by the hierarchy 



m t ^> m t v (three-momentum, "soft" scale) 3> m t v 2 (kinetic energy, "ultrasoft" scale) 3> Aqcd 



(1) 



The particle-antiparticle propagation is described by the terms in the effective theory Lagrangian which are bilinear 
in the top quark and antitop quark fields, 



C(x) 



Sm t >ip P (x) + (ip P (x) -> Xv( x )) ■ 



(2) 



where the fields i/j p and \p destroy top and antitop quarks with soft three-momentum p in the ti cm. frame and T t 
is the on-shell top quark decay width. The term Smt is a residual mass term specific to the top quark mass definition 
that is being used; for our analysis we employ the IS mass scheme 0,0]. 

Up to NLL order the top-antitop quark pair interacts only through the effective Coulomb potential [lsf . 



V c (p,q) 



A-TxCFOL a (jritv) 



a s (m t v) 

47T 



-A) In 



a i 



(3) 



where k = p — q is the momentum transfer and /3q = 11/3(7^ — A/STrif is the one-loop QCD beta function, 
di = 31/9Ca — 20 /9Tnf the coefficient of the one- loop correction to the effective Coulomb potential, and Ca = 
2> 7 Cf = 4/3, T =1/2 are SU(3) group theoretical factors. For the number of light flavors we take rif = 5. The 
parameter v is the vNRQCD renormalization scaling parameter used to describe the correlated running of soft and 
ultrasoft effects in the effective theory. Here, v = 1 corresponds to the hard scale at which the effective theory is 
matched to the full theory, and v = vq, vq being of the order of the typical ti relative velocity, is the scale where 
the matrix elements arc computed. The evolution of the Wilson coefficients from the matching scale down to the 
low-energy scale sums logarithms of the velocity to all orders and is governed by the velocity renormalization group 
equations [lo| . 
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Top-antitop quark production in the nonrelativistic regime in the LL and NLL approximation in a 3 Si spin triplet 
or a ! 5o spin singlet state is described by the currents 



J lp = ippO'ji.i^X- 



(4) 



where c\.j{v) and Co{y) are the corresponding Wilson coefficients. The currents do not run at LL order but UV- 

lead to non-trivial anomalous dimensions at NLL order 



divergences in effective theory two-loop vertex diagra ms id 

BET 



which result in a scaling of the Wilson coefficients 13, 

djM=cij(l) exp(/(^S 2 =2)) , cq{v) =co(1) exp(/(^,S 2 - 0)) , (5) 

where S 2 is the square of the total ti spin. We have adopted the convention that the matching conditions at v = 1 
only account for QCD effects, so at LL order we have Ci(l) = cq(1) — 1. The NLL order QCD matching conditions 
relevant for e + e~ — > tiH in the large Higgs energy endpoint region are discussed in Sec. OH 

Through the optical theorem the ti production rate for a ti invariant mass Q 2 ~ Am 2 involves the imaginary part 
of the time-ordered product of the production and annihilation currents defined in Eqs. 



A l f{Q\m u v) 



Aq(Q 2 ,m u v) = i 



E 

p,p' 

E 

pp' 



T J[] p ,(0)Ji P (x) 0) = 2N c 6 lk G c (a,v,m t ,v), 



d x i 



TJ T p ,(0)Jo,p(x) 0) = N c G c (a,v,m u v), 



(6) 



(7) 



where v — ((y^Q 2 — 2m t — 28mt + iTt)/m t )^ is the cm. top quark (effective) relative velocity and q = (^/Q 2 — 2m t , 0). 
The term G c is the zero-distance S-wave Coulomb Green function of the nonrelativistic Schrodinger equation with 
the potential in Eq. J3J . To compute the Green function we use the numerical techniques and codes of the TOPPIC 
program developed in Ref. [2ol | and determine an exact solution of the full NLL Schrodinger equation following the 
approach of Refs. 01 • 



3. Comments on the Computation 

The LL vNRQCD result for the Higgs energy distribution, including the QCD effects coming from the Coulomb 
potential in Eq. (J3J, the finite top quark lifetime, and the Wilson coefficients of the currents in Eqs. (JSJ, is given by 

-^&B*Eir) - SNC[{1 + XH ^~ AXH]1/2 {^)Fi + c 2 MFT Z ) InrlO,^,,)] . (8) 

For details on the explicit calculation and the definition of the formfactors Fj see |9(- We note that the NLL (0(a s )) 
matching conditions for the three triplet Wilson coefficients C\j depend on the ti spin configuration (i.e. on j) since 
the kinematic situation for tiH production in the large Higgs energy endpoint is not invariant under separate rotations 
of the spin quantization axis. However, for our purposes it is sufficient to define a triplet Wilson coefficient that 
is averaged over the three spin configurations. Using such an averaged triplet Wilson coefficient the Higgs energy 
spectrum at NLL order can also be cast in the simple form of Eq. ijHJ . 

At NLL order, we need to account for the 0(a 2 ) contributions to the Coulomb potential in Eq. the NLL 
running of the coefficients c\ and cq, and their 0(a s ) matching conditions at v = 1. The latter hard QCD corrections 
are process specific and cannot be inferred from results obtained in earlier computations for the ti threshold in 
e + e~ — * ti. We have extracted the matching conditions from the codes for the Standard Model amplitude for 
e + e _ — > tiH provided in Ref. 0. With the ansatz 

Cqa{v = l,v r s,m t ,m H ) = 1 + < ~' F ° Ls — — ^c ,i(\/s, mt, msr) (9) 

we have determined the 0(a s ) matching conditions numerically by matching the 0(a s ) vNRQCD prediction at 
fi = m t {v = 1) to the full theory results close to the large Higgs energy endpoint. The relative uncertainties for this 
numerical procedure are below 1% for <5co,i. 
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Figure 2: Diagrams describing virtual one-loop QCD corrections in the Standard Model to compute the NLL matching 
conditions for the operators Ji,p and Jo,p that describe ti production for invariant masses Q 2 « 4m?. Self energy diagrams 
are implied. 
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Figure 3: Higgs energy spectrum at LL (dotted lines) and NLL (solid lines) order in the nonrelativistic expansion for the 
vNRQCD renormalization parameters v = 0.1, 0.2, 0.4 and at the Born level and at 0(a s ) for ll = y/s, y/s\v\ for the parameters 
(a) yfe = 500 GeV, m H = 120 GeV and (b) yfl = 550 GeV, m H = 140 GeV. The top IS mass has been set to mf = 180 GeV 
and the other parameters are Vt = 1.55 GeV, Mz = 91.1876 GeV, M w = 80.423 GeV, c w = M w /M z and a' 1 = 137.034. 



4. Numerical Analysis 

In Figs.|21the predictions of the Higgs energy spectrum in the full kinematic range are displayed for two cases. In 
the large Higgs energy endpoint region we have shown the LL (dashed lines) and NLL (solid lines) results in the 
nonrelativistic expansion. (See the figure captions for details on choice of parameters.) At LL order the upper (lower) 
curve corresponds to v = 0.1 (0.4), while at NLL order the upper (lower) curve corresponds to v = 0.2 (0.1). The 
curves show the typical behavior of the prediction of the nonrelativistic expansion for any choice of parameters. While 
the LL predictions have a quite large renormalization parameter dependence at the level of several tens of percent, 
the NLL results are stable. Here, the variation due to change of the renormalization parameter is around 5%. The 
stabilization with respect to renormalization parameter variations at NLL order arises mainly from the inclusion of 
the 0(a s ) QCD corrections to the Coulomb potential, Eq. J2J). Moreover, the NLL curves lie considerably lower than 
the LL ones. This behavior is well known from the predictions for e + e _ — > ti at threshold and originates from 

the structure of the large negative 0(a s ) QCD corrections to the Coulomb potential, Eq. ©, and from the sizeable 
negative QCD corrections to the matching conditions in Eq. 10. 

In principle this behavior is a point of concern because it could indicate that the renormalization parameter 
variation might be an inadequate method to estimate theoretical uncertainties. Fortunately, the top quark mass is 
sufficiently large such that the regions where the conventional fixed-order expansion (in powers of the strong coupling) 
and where the nonrelativistic expansion (described by the effective theory) can be applied are expected to overlap. 
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[GeV] 


m H [GeV] 


cr(Born) [fb] 


ff(a.) [fb] 


cr(NLL) [fb] 


cr(NLL) 
CT (Bo„) 


a (NLL) 


(7 (NLL) |^ | <o.2 


500 


120 


0.151 


0.263 


0.357(20) 


2.362 


1.359 


1.78 


550 


120 


0.984 


1.251 


1.342(37) 


1.364 


1.073 


1.66 


550 


160 


0.134 


0.207 


0.254(12) 


1.902 


1.226 


1.74 


600 


120 


1.691 


1.939 


2.005(30) 


1.185 


1.034 


1.66 


600 


160 


0.565 


0.700 


0.745(18) 


1.319 


1.065 


1.68 


700 


120 


2.348 


2.454 


2.485(13) 


1.058 


1.012 


1.68 


700 


160 


1.210 


1.303 


1.328(11) 


1.098 


1.020 


1.69 



Table I: Collection of cross sections and K factors for various cm. energies and Higgs masses and top quark mass m t = 
180 GeV. Here, a(Born) refers to the Born cross section and a(a s ) to the 0(a s ) cross section in fixed-order perturbation 
theory using fi = \/s as the renormalization scale, the choice employed in the analysis of Ref. [(J. The term a (NLL) refers to 
the sum of the 0(a s ) fixed-order cross section for v > 0.2 using fi = %/sv and the NLL nonrelativistic cross section for \v\ < 0.2 
with the renormalization parameter v — 0.2. (See Figs.[3]where the vertical lines indicate the Higgs energy where v — 0.2.) 

To demonstrate this issue we have also displayed in Figs.[3] the fixed-order (i.e. without summation of Coulomb 
singularities) predictions at the Born (dashed lines) Q and the 0(a s ) level (solid lines) The two 0{a s ) curves 
correspond to the renormalization scales fi — y/s (lower curves) and /i = \y/sv\ (upper curves), where v is the tt 
relative velocity defined below Eq. 10. The latter choice for the fixed-order renormalization scale is motivated by the 
fact that the relative momentum of the top pair is the scale governing the Coulomb singularities contained in the 
fixed-order expansion close to the large Higgs energy endpoint. This choice for the fixed-order renormalization scale 
is therefore the more appropriate one near the Higgs energy endpoint. The results in Figs. [3] demonstrate the overlap 
between the 0(a s ) fixed-order prediction and the NLL nonrelativistic one in the region where the it relative velocity 
is approximately 0.2. (The Higgs energy with v = 0.2 is indicated in each panel by the solid vertical line.) The overlap 
improves for increasing cm. energies or decreasing Higgs masses. This indicates that in the overlap regions the higher 
order contributions summed in the nonrelativistic prediction and the higher order relativistic corrections contained 
in the fixed-order result are both small. For smaller cm. energies or increasing Higgs masses, on the other hand, the 
NLL nonrelativistic predictions tend to lie slightly above the 0(a s ) fixed-order results (for fi — y/sv) illustrating the 
impact of the higher order corrections to each type of expansion. The discrepancy, however, remains comparable to 
the uncertainties estimated from the renormalization parameter variation of the NLL nonrelativistic prediction. We 
therefore conclude that the renormalization parameter variation of the NLL order nonrelativistic prediction should 
provide a realistic estimate of the theoretical uncertainties in the large Higgs energy region. The results in Figs.EH 
also demonstrate that the region of parameter space where the top quark pair is nonrelativistic increases for smaller 
cm. energies (or larger Higgs masses). 

Let us now discuss the numerical impact of the nonrelativistic contributions in the large Higgs energy region on the 
total cross section. In Tab.[|]the importance of the summation of the Coulomb singularities and the logarithms of the 
top quark velocity is analyzed numerically for various choices of the cm. energy and the Higgs mass. For all cases the 
top quark mass mj s = 180 GeV is used and the other parameters are fixed as in Figs. 03 See the caption for details 
on the various entries. The numbers for cr(NLL), which are determined from combining the NLL nonrelativistic 
predictions in the Higgs energy end point for \v\ < 0.2 with the fixed-order 0(a s ) prediction for smaller Higgs 
energies with \v\ > 0.2, represent the currently most complete predictions for the total cross section of the process 
e + e~ — > tiH as far as QCD corrections are concerned. For cr(NLL) we have also given our estimate for the theoretical 
error. For the fixed-order contribution (v > 0.2) we have estimated the uncertainty by taking the maximum of the 
shifts obtained from varying fj, in the ranges [y/s, 2^/s], [y/s, y/s/2] and [y/sv, y/s\; for the nonrelativistic contribution 
in the end point we have assumed an uncertainty of 5% for all cases. For the numbers displayed in Tab. [i] both 
uncertainties were added linearly. 

The results show that the enhancement of the cross section due to the summations in the large Higgs energy 
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region is particularly important for smaller cm. energies and larger Higgs masses, when the portion of the phase 
space where the nonrelativistic expansion has to be applied is large. Here, the higher order summations contained in 
the nonrelativistic expansion can be comparable to the already sizeable 0(a s ) fixed-order corrections and enhance 
the cross section further. This is advantageous for top Yukawa coupling measurements for the lower cm. energies 
accessible in the first phase of the ILC experiment. For higher cm. energies the effect of the nonrelativistic summations 
of contributions from beyond 0(a s ) is less pronounced and decreases to the one-percent level for cm. energies above 
700 GeV. For all cases, except for very large cm. energies around 1000 GeV, however, the shift caused by the terms 
that are summed up in the nonrelativistic expansion exceeds the theoretical error 

In Tabled in the last column, the ratio of the NLL nonrelativistic cross section and the 0(a s ) fixed-order cross 
section (with the approximation rnj? ole = ml s ) for |u| < 0.2 is also shown. Interestingly, the higher order summations 
lead to correction factors ranging between about 1.7 and 1.8 that are only very weakly dependent of the cm. energy 
and the Higgs mass. This fact might prove useful for rough implementations of nonrelativistic tt effects in other high 
energy processes. 
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